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A3STRACT 


Expressions  for  a  generalized  Modified 
Transitional  3utter worth-Chebyshev  (MTBC)  filter  are 
derived.  The  characteristics  of  this  filter  as  applied 
to  digital  filter  design  are  investigated.  It  is  shown 
that  by  adjusting  location  and  order  of  the  inserted 
zeros,  the  cut-off  slope  rate  of  the  filter  can  be 
traded  for  maximum  attenuation  in  the  stop-band. 

The  performance  of  this  MTBC  filter  is  compared  to 
that  of  Butterworth,  Chebyshev,  transitional 
Butterworth-Chebyshev  filters  together  with  those 
suggested  by  other  investigators  [1]-[3].  It  is  shown 
that  the  step-band  attenuation  can  be  significantly 
increased  without  great  sacrifice  of  cut-off  slope 
rate. 

Step  response  of  this  MTBC  filter  is  also  obtained 
and  compared  with  other  filters.  Various  tabulations 
as  well  as  graphs  of  this  filter  are  given  for  design 
purposes.  A  computer  program  is  developed  for  the 
design  of  this  filter. 
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I.   INTRODUCTION 


The  digital  filter  is,  as  defined  by  Rabiner  et  al  [10], 
"  a  computational  process  or  algorithm  by  which  a  digital 
signal  or  sequence  of  numbers  (acting  as  input)  is 
transformed  into  a  second  sequence  of  numbers  termed  the 
output  digital  signal". 

The  area  of  digital  filtering  can  be  divided  into  two 
major  subdivisions  as  Finite  Impulse  Response  (FIR)  filters 
and  Infinite  Impulse  Response  (IIR)  filters. 

During  the  development  of  digital  signal  processing,  the 
interest  cf  the  investigators  in  IIR  and  FIR  filters  varied. 
Before  the  introduction  of  the  FFT  algorithm  by  Cooley  and 
Tukey  (1965)  IIR  filters  were  much  more  efficient  than  FIR 
filters.  Stccham's  work  [13]  on  the  FFT  method  of 
performing  convolution  indicated  that  implementation  of 
high-order  FIR  filters  could  be  made  extremely 
computationally  efficient  ;  thus,  comparison  between  FIR  and 
IIR  filters  are  no  lcnger  strongly  biased  toward  the  latter 
[5].  Because  FIR  filters  require  very  high  orders  to 
produce  a  sharp  attenuation  shape,  they  are  not  often  used 
for  real-time  filtering  of  waveforms.  Recently,  due  to  the 
increase  in  computing  capabilities  in  digital  signal 
processing  and  the  availibity  of  long  charge  transfer  device 
(CTD)  tapped  delay  lines  (TDL) ,  FIR  filters  are  favored  over 
IIR  filters.  However,  in  applications  like  design  of 
digital  ccmb  filters  IIR  filters  are  the  unique  alternative. 

There  are  three  basic  design  techniques  of  IIR  digital 
filters  [5]. 
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First  method  is  the  direct  design,  which  is, 
appropriately  placing  poles  and  zeros  to  approximate 
required  frequency  response. 

A  second  method  is  to  use  an  optimization  procedure  to 
place  the  foles  and  zeros  to  match  arbitrary  frequency 
response  specifications. 

Finally  the  third  technique  makes  use  of  highly  advanced 
art  of  continuous  filter  design.  This  technique  of 
designing  digital  filters  from  continuous  filters  by  means 
of  mathematical  transformations  is  the  most  popular  IIR 
digital  filter  design  technique. 

Standard  Z-transform, Bilinear  Z-transform  and  the 
matched  Z-transform  make  possible  direct  transformation  from 
S-domain  to  Z-domain,  preservinq  essential  characteristics 
of  analog  frequency  response. 

Existence  of  frequency  transformations  reduces  the 
problem  to  design  a  frequency  normalized  prototype  low-pass 
filter. Then  using  appropriate  frequency  transformation,  this 
prototype  may  be  converted  into  desired  band-pass, 
band-reject  cr  high-pass  filter.  Popular  prototype  filters 
are  Butterworth,  Chebyshev,  elliptic  and  hybrid  transitional 
filters.  Frequency  transformations  for  digital  filters  are 
discussed  in  various  literatures  (  [6]  and  [7]  ). 

The  problem  of  designing  low-pass  prototype  filters, 
which  possesses  better  stop-band  attenuation  and  cut-off 
slope  characteristics  than  existing  prototypes  has  always 
attracted  the  researchers  in  the  signal  processing  area. 

3udak  and  Aronhime  suggested  [1]  modification  of 
maximally  flat  rational  functions  by  introducing  a  pair  of 
finite  transmission   zeros   such   that   the   maximally   flat 
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characteristic  is  maintained  but  the  cut-off  slope  can  be 
made  steeper  without  great  sacrifice  of  stop-band 
attenuation. 

Dutto  Roy  [2]  investigated  a  more  general  case  allowing 
insertion  of  multiple  pairs  of  transmission  zeros,  either 
coincident  or  distinct. 

Introducing  multiple  pairs  of  jw-axis  zeros  in  all  pole 
Chebyshev  transfer  functions  are  investigated  by  Agarwal  and 
Sedra  [  3  ]. 

The  most  attractive  feature  of  these  finite  zero  filters 
is  that  they  offer  the  filter  designer  a  great  degree  of 
freedom  in  choosing  the  location  and  order  of  the  zeros  to 
trade  cut-off  slope  for  stop-band  attenuation. 

In  this  thesis,  a  modified  Transitional 
Butterworth-Chebyshev  filter  is  developed,  which  is  a  more 
general  case,  introducing  finite  coincident  or  distinct 
multiple  pairs  of  transmission  zeros  in  transitional 
Butterworth-Chebyshev  filter. 

Trade-off's  between  tne  order  of  the  filter,  the  order 
of  transmission  zeros,  stop-band  attenuation  and  cut-off 
slope  are  pointed  out.  Graphs  helpful  in  the  design  of  such 
filters  are  obtained. 

Performances  of  Butterworth,  Chebyshev,  Transitional 
Butterworth-Chebyshev  filters  and  the  designs  suggested  in 
references  [2]  and  [3]  are  compared  with  those  represented 
in  this  thesis  for  the  orders  of  three  through  eleven.  A 
computer  program  is  developed  to  implement  the  filters 
mentioned  above. 

In  addition,  the  time-domain  response  of  digital  filters 
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is  studied.  There  are  many  applications,  such  as  digital 
HTI  filters,  for  which  one  is  interested  in  the  transient 
responses  cf  filters  that  are  specified  in  the  frequency 
domain.  Step  responses  of  the  filters  that  are  discussed  in 
this  report  are  plotted  and  compared. 
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II.   DERIVATION  OF  MODIFIED  TRANSITIONAL 
BOTTERWORTH-CHEBYSHEV  FILTERS 


A.   INTRODUCTION 


The  most  popular  technique  for  designing  IIR  digital 
filters  is  tc  digitize  an  analog  filter  that  satisfies  the 
design  specifications  [5].  There  are  many  techniques  for 
designing  analog  low-pass  prototype  filters.  Among  the  well 
known  analog  filter  classes  are  the  maximally  flat 
(Butterworth)  and  equal  ripple  (Chebyshev)  filters. 

Butterworth  filters  are  simple,  excellent  in  the 
pass-band  and  monotonic  in  both  pass-band  and  stop-band.  The 
Chebyshev  filters  are  superior  at  and  near  cut-off  frequency 
and  at  stcp-cand. 

The  transitional  Butterworth-Chebyshev  (TBC)  filters 
combine  the  desirable  attributes  of  these  two  filters  in  a 
single  approximation  that  is  given  by 

where 

k=  Weighting  factor 

n=  Order  of  filter 

C  .  (w)  =  (n-k)—  order  chebyshev  polynomial. 
When   k  =  n,   |F(jw)|2   is   identical   with   the   Butterworth 
function.  When  k=0,  JF(jw)|2  is  identical  with  the  Chetyshev 
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function.  With  a  varying  value  of  k,  a  TBC  filter  possesses 
seme  characteristics  of  each.  As  k  approaches  to  n,  a  TBC 
filter  behaves  more  like  a  Butterworth  filter,  as  k 
approaches  to  zero,  it  behaves  more  like  a  Chebyshev  filter. 

In  this  chapter,  modification  of  TBC  filters  by 
introducing  finite  coincident  or  distinct  multiple  pairs  of 
transmission  zeros  will  be  discussed.  It  will  be  shown 
that,  using  a  weighting  factor  k,  and  the  location  and  order 
of  inserted  zeros  as  parameters,  attenuation  in  the 
stop-band  may  be  traded  for  sharpness  of  the  cut-off 
characteristics. 

Expressions    for   the  cut-off   slope   and   minimum 

attenuation  in  the  stop-band  are  derived  in  terms   of   order 

of  the  filter,  weighting  factor,  location,  and  order  of 
inserted  zeros. 


B.   MODIFICATION  OF  THE  TBC  FILTERS  WITH  COINCIDENT 
TRANSMISSION  ZEROS 


Introducing   m   identical  pairs  of  transmission  zeros  at 
±j«0  to  eg.  (2-1)  ,  we  have 


In  order  to  normalize  |F(jw)|*,  i.e.  to  force  |F(jw)|*  to  be 
egual  to  1/2  at  m= 1 ,  the  constant  K  should  be 
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Then    eq.     (2-1)    becomes 


F<j«) 


2m. 


C*S-Mix)      +  <W-t)2,«  ^  cn\kM 


w 


where  n>2m,  because  of  the  low-pass  characteristics   of   the 
function,  and  the  Chebyshev  polynomial  C^  (w)  is  defined  by 


Cos  (n  Cos"' to) 


,  toi  I 


c>)  = 


(«) 


2r\-l      In. 


eos^Cn&slr'w)^  ***«?*   ,  w)i   (J-6) 


1 •       Slope   at    cut-off   frequency 


Substitutinq    eq.     (2-5)    into   eq.     (2-4),    we   obtain 


FCjw) 


(Wo*-  U3r) 


2m 


Let 


h(co)  =     u^2t  ^[("-It)  Cos'io] 

Usinq  these  values,  eq .  (2-7)  becomes 


■ffa) 


F(>)l  "  {(«,)  + (w,M)*m  w(«) 


(a-.o) 
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Taking  the  derivative  of  eg.  (2-10)  with  respect   to   w,   we 
get 


M 


At   the    cut-off    frequency,    w=1,    we    have 

f(w)  =  (*.*- i)2"" 

kfo)-  1 


(w0 


Substituting   these  values  into  eg.  (2-11)  and  simplif iying, 
we  obtain 


f(y*)\    - 


m 


JT  (uo^-O 


+ 


fe.4  (n-lp 


H 


For  k=n  ,  the  result  agrees  with  the  cut-off  slope  of  MB 
function,  which  is  derived  in  [2].  For  k=0,  the  result 
agrees  with  the  cut-off  slope  of  MC  function  as  derived  in 
[3]. 

2 •   St o^^band  characteristics 


In  general,  stop-band  characteristics  of  finite  zero 
filters  will  be  of  the  form  shown  in  FIG.1*. 


*  Eecause  of  the  inherent  limitations  in  the  plotting 
subroutine  utilized  to  provide  the  graphs  in  this  thesis,  it 
was  necessary  to  add  supplementary  axes  to  show  proper 
scaling  for  seme  of  the  graphs. 
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Figure    1    - 


GENERAL  STOP-BAND  CHARACTERISTICS  OF  THE  FINITE 
ZERO  FILTERS 
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The  zeros  introduced  will  cause  a  peak  in  the  stop-*band,  at 
a  freguency  vy>w0  .  The  minumum  attenuation  in  the  stop-band 
may  be  defined  as 


c<*T8c  -  JO  log  Tp 


(hi) 


where 

i 

i 

Ff- 

FCju)f) 

In    th 

e   stop-t 

and   eg. 

(2-4)     becomes 

H 


t^) 


(u)ol-u)/) 


2»n 


(u30*.^)2^^  (Oo1--!)2"  S^-1    ^ 


(2-20) 


Let 


f  l*r)  = 

Vi  (u)f)  = 


2n\ 


(uJo'-U).1) 


(«f) 

(2-22) 


Taking   derivatives   of   these   values  and  substituting  them 
into  eq.  (2-11)  results  in 


F(j^) 


f  =  V  n-^ 


Thus 


1  60^  )2m  ^iU"kVz  ^  (^^  [- w-^w)  *»  ] 

(^-23) 

(i-2v) 


Combining  eg.  (2-19)  and  eg.  (2-23),  we  obtain 


u), 


f,= 


1+  (u0o-<;       o2 


a' 


,      „     Nf\-2m  2.rr\ 


U) 


l(r\-2m) 


H 
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And  minunium  stop-hand  attenuation  will  be  given  by 

CX^tsc   -     10  n  log  ( H \  +i0ry\   I09  (u)0) 


+,20  m  io 


3 


^  A-2rr\  > 

I            -p     BW    V            ., 

n-£m 

10<M 

n 

U3oz 

f  6(n-*v\-te->t)    6-^ 


Plots  of  stop-band  attenuation  and  cut-off  slope  of 
MTEC  filters  with  two  coincident  transmission  zeros,  for 
orders  3  thrcugh  1 1  are  given  in  figures  2  and  3. 


C.   MODIFICATION  OF  T3C  FILTEHS  WITH  DISTINCT  TRANSMISSION 
ZZBCS 


Consider   the   n  —  order   TBC   function   with  in  pairs  of 
finite  distinct  zeros  at±jw«  ,  where  1=1,2,...  ,1a  and  w-  >1. 

The   rragnitude   squared   function   of   TBC  function  with 
distinct  transmission  zeros  will  be  given  by 
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Kji) 
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7f  ( u>-.1- u>*) l  +  J  c «*- o    ^2*  c0.k w 
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•  =  i 


Let 


•vi 


IS  I 

Putting  these  values  into  eg.  (2-27)  results  in 
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F(jw) 
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Figure  2  - 


ST0P-3AND  ATTENUATION  OF  MTBC  FILTER  WITH  TWO 
COINCIDENT  ZEROS 
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Figure    3    -       COT-OFF    SLOPS    OF    MTBC    FILTER    WITH    TWO    COINCIDENT 

ZEROS 
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1 •   Sto£-band  and  cut-off  characteristics 


The  frequency  wp   at  which  the  stop-band   peak   will 
occur,  is  given  by 


FCj*) 


=  0 


H 


Ui=rU), 


Using  eq.  (2-6)  and  taking  derivative   of   eg. 
respect  tc  w,  we  obtain 


(2-28)       with 


KjuO 


JL\H$>)\  Wio)[u  560)  a 


K(ta) 


1 


(«.) 


Combining  eg.  (2-29)  and   eg.   (2-30)   the   equation   to   be 
solved  fcr  w_   may  be  found  to  be 


a  10 


f 


m 


i 


(i0;z-u)pz) 


n  =  0 


(«i) 


which  agrees  with  eg.  (2-24)  when  w;  »s  are  identical.  Using 
real  solutions  of  eg.  (2-31),  stop-band  attenuation  peaks 
may  be  fcund  to  be 


£<-.  -  JO  iog  tjo 


ft: 


where 


Fp;  =  |F(jioPl-) 


and  minimum  stop-band  atttenuation  is  given  by 


c< 


mm.  =  RtaT* 
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(is*) 
(4*3) 


(a.*y) 
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Using   eg.  (2-5)  and  differentiating  eg. (2-27)   with  respect 
to  w,  the  cut-off  slope  becomes 


FQ«)     — 


I 


ZU? 


k  +  Ck-^)2  +  £ 


m 


«=i 


10;*-! 


H 


which  agrees  with  eg.  (2-17)  when  the  w-,  »s  are  identical. 

Plots  of  cut-off  slope  and  stop-band  attenuation  of 
MTBC  filters  with  two  distinct  transmission  zeros,  for 
orders  3  through  11  are  given  in  figures  3-17. 


D.   SUMMARY 


The  Transitional  3utterworth-Chebyshev  filter  combine 
the  best  features  of  the  Butterworth  and  Chebyshev  filters. 
A  Modified  Transitional  Butterworth-Chebyshev  filter, 
obtained  by  introducing  finite  transmission  zeros  to  a 
Transitional  Butterworth-Chebyshev  filter,  possesses  cut-off 
slope  and  stop-band  attenuation,  which  are  dependent  on  the 
modification  parameters  w0  and  m.  The  closer  the  wa  is  to 
unity,  the  steeper  the  cut-off  slope.  However,  this 
improvement  in  the  cut-off  region  results  in  degradation  of 
stop-band  attenuation.  Thus,  using  w0  and  m  as  parameters, 
an  advantageous  trade  between  attenuation  in  stop-band  and 
sharpness  of  the  cut-off  characteristics  can  be  made. 
Graphs  are  given  to  serve  as  guides  in  trading  cut-off  slope 
for  stop-band  attenuation. 
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Figure  a  -   CUT-OFF  SLOPE  OF  MTBC  FILTER  WITH  TWO  DISTINCT 

ZEROS  <w=1«06,  n=5> 
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Figure    7    -      CUT-OFF    SLOPE    OF    MTB-C    FILTER    WITH    TWO    DISTINCT 

ZEROS    <w    =1«06>,    n=8> 
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Figure    8    -      CUT-OFF    SLOPS    OF    MTBC    FILTER    WITH    TWO    DISTINCT 

ZEROS    <w    =1«06,    n=9> 
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Figure  9  -   CUT-OFF  SLOPE  OF  HTBC  FILTER  WITH  DISTINCT 
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Figure  10  -   CUT-OFF  SLOPE  OF  *TBC  FILTER  WITH  DISTINCT 
TRANSMISSION  ZEROS  <w  =1*06  ,n-11> 
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Figure  11  -   STOP-BAND  ATTENUATION  OF  MTBC  FILTERS  WITH  TWO 
DISTINCT  TRANSMISSION  ZEROS  <w  =1*06,  n=5> 


34 


1.3 


1.0 


i-a 


k-o 


r.o 


U> 


6 


cO 

d 

o 

-r 


-5 

c 

1      : 
Q 


C] 


U3a  (rad/s<2cj  - 

K-SHRLF.    1-  CVr-0L   UN 


-3C  D|_p-„  i„  GGE-'-Oi 


!ki!Ts 


r  i  ; r- i_ 


!  .  s_i!  •  » 


\ 

.  \  \ 

'      \ 

\N   N 

\\\  \ 

\  \ 
\\  \ 

■\\\\    \ 

\\V\  s 


%\\ 


Figure  12  -   STOP-BAND  ATTENUATION  OF  MTBC  FILTERS  WITH  TWO 
DISTINCT  TRANSMISSION  ZEROS  <w  =1*06,  n=6 > 
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Figure  14  -   STOP-BAND  ATTENUATION  OF  I1TBC  FILTERS  WITH  TWO 
DISTINCT  TRANSMISSION  ZEROS  <v  =1«06,  n=8> 


37 


*o 


1.0 


3-0 


H»o 


r.o 


Co 

i 


5 

i 


I   I 


C0 

o 

e 

5  ' 

c 
<i 

J> .. 

I   ? 

4 

<0 


0)i  (racL  jse.c) 


U),  =  <.o6 
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Figure  16  -   STOP-BAND  ATTENUATION  OF  MTBC  FILTERS  WITH  TWO 
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Figure    17    -      STOP-BAND    ATTENUATION    OF    MTBC    FILTERS    WITH    TWO 
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III.   COMPARISON  OF  MrBC  FUNCTION  WITH  B,tMBxCJ.MC  AND  TBC 

FUNCTIONS 


A.   INTRODUCTION 


In  chapter  2,  derivation  of  MTBC  filter_is  given.  It  is 
shown  that,  stop-band  attenuation  and  cut-off  slope  rate  of 
this  filter  depend  on  modification  parameters  m,  k,  and  w0  . 
So  the  filter  designer  have  the  flexibility  to  trade  cut-off 
slope  for  stop-band  attenuation  by  changing  these 
parameters,  without  changing  the  order  of  the  filter. 

In  this  chapter,  performance  of  MTBC  filter  is  compared 
with  B,  MB,  C,  MC,  and  TBC  filters. 

Formulas  for  stop-band  attenuations  and  cut-off  slopes 
of  all  these  filters  are  given  in  Table  I. 


B.   MTBC  FUNCTION  V.S.  B  FUNCTION 


Cut-off  slope  of  MTBC  function  is  given  by 
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From  eg.  (3-1)   ratio  of   cut-off   slopes   of   these   two 
functions  may  be  found  as 


'MT6C./8   -    "? 
06 


=■  fl-V- 


J2 


m 


n(\0oM) 


(••0 


rMi8c/&  v.s.  w0  for  n=3  through  11  and  ra=  1,2  are  given  in 
Figures  18  and  19.  Stop-band  attenuation  of  MTBC  may  be 
written  as 


^mtgc  =   °^q  +  6(n-m-fe-i)  ■+  o20  mr\  Loo 


n-im     u30x-4 


n 


lOo2" 


(3-0 


Then   the  difference  between  stop-band  attenuations  is  given 
by 


dMT8c/s=  <**t«c  -  °^g   -  £Gi-m-k-l)+AOmloq 


A-2yy\       lQoa--l 

n         too1 


H 


Plots      of      d 


HTtC 


/Q   v.s.  w0   for  n=3  through  11  and  m=1f2 


are  given  in  FIG  20  and  21. 

Figures  19-21  indicate  that,  MTBC  filter  can  be  made  10 
times  steeper  than  B  filter  still  having  40  dB  more 
attenuation  at  the  stop-band. 


C.   MTBC  V.S.  MB  FUNCTION 


Using  Table  I,  cut-  off  slope   ratios   of   MTBC   and   MB 
filters  may  be  found  as 


/         5m  +n  (lOoM; 


(«) 


Plcts   of    ^HliCJMa    v.s.    w0      for   n  =  3    through    11    and      m=1,2      are 

given    in    figures   22    and   23. 

The    difference   between   the    stop-band      attenuations      of      MTBC 
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and  MB  is  given  by 


Equation  (3-7)  indicates  that,  stop-band  attenuation 
difference  between  these  filters  doesn't  depend  on  w0 .  For 
n=10,  i=1,  k=1,  with  48  dB  more  attenuation,  cut-off  slope 
ratio  may  be  changed  from  2.5  to  8.0  by  changing  w©  from 
1.06  to  2.4. 


D.   MTBC  FUNCTION  V.S.  C  FUNCTION 


The   ratio   of  cut-off  slope  of  MTBC  function  to  cut-off 
slope  of  C  function  is  given  by 


MTBC 


:/c  =■  i  + 


im 


k(l*+l-i*) 


A2(ul»l-0 


rV 


Plot  of  rMTg<i/c  v.s.  w„  for  n=3  through  11  and  m=1  ,2  are 
given  in  Figures  24  and  25.  The  difference  between 
step-band   attenuations   of   MTBC   and   C   functions   is    given   by 

- 

Plots  of  &mtcJc  v*s  w°  ^or  n=~*  "through  11  and  m=1,2  are 
given  in  Figures  26  and  27.  The  Chebyshev  filter  is  known 
to  provide  much  steeper  cut-off  slope  than  the  corresponding 
B  and  TBC  filters.  Figure  25  shows  that  the  cut-off  slope 
of  the  MTBC  filter  can  be  made  1.4  times  steeper  than  that 
of  the  Chebyshev  filter,  with  n=5,  m=2,  and  w„=1.Q6. 


E.   MTBC  FUNCTION  V.S.  MC  FUNCTION 
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The    ratio   of   cut-off   slope   of      MTBC      filter     to      cut-off 
slope   of    MC    filter    is    given    by 


(3..0) 


Plots   of   Zhiqc/mc     v,s*    w«      for   n=3    through    11    and      m  =  1,2      are 
given    in    Figures   28   and  29. 

The    difference    between   stop-band  attenuations  of   aTBC    and   MC 
functions   is    given    by 

For    a   given    k,    ^UTar/MC     "*"s    a^-^ways   constant. 


F.   MTBC  V.S.  TBC  FUNCTIONS 


The   ratio   of   the  cut-off   slopes   of   MTBC   and   TBC 
functions  is  given  by 


^T6c/tbc  ■!  + 


m 


H 


0o.l-J)  [n2-t  k(k-v  l-4<0] 

Plots   of  ^tcncljzc   v«s»  w»   f°r  n=3  through  11  and  m  =  1  ,2   are 
given  in  Figures  30  and  31. 


The   difference   between   the   stop-band  attenuations  of 
these  two  functions  is  given  by 


^M-rac/rac  =  -  6m  -*■  3.0  w   Ua 


n-£rvi        lo*M 


r\ 


Wo1 


(«-b) 


Plots      of      3MTe£/Tec  v's«    wo       for   n=3    through    11    and   m=1,2 
are    given    in    Figures   32   and    33. 
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G.   SUMMARY 


Using  the  location  and  the  order  of  the  inserted  zeros 
and  the  weighting  factor  as  parameters,  the  characteristic 
curve  of  the  MTBC  filter  can  be  made  steeper  than  that  of 
the  conventional  all-pole  filters  without  greatly 
sacrificing  either  stop-band  attenuation  or  flatness  in  the 
pass  band.  A  Modified  Chebyshev  filter  provides  slightly 
better  performance  than  the  MTBC  filter,  at  the  expense  of 
sutstantial  degradation  of  the  pass-band  flatness.  Graphs 
are  given  to  help  in  the  comparison  and  in  determining  the 
numerical  advantages  gained  by  increased  complexity. 
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RATIO    OF    CUT-OFF    SLOPES    OF    MTBC    AND    B    FUNCTIONS 
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Figure     19    -      RATIO    OF    CUT-OFF    SLOPES    OF    MTBC    AND    B    FUNCTIONS 

V«S«    W    <M=2> 
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Figure    20    -       THE    DIFFERENCE    BETWEEN    THE    STOP-BAND 
ATTENUATIONS    OF    21TBC    AND    B    FUNCTIONS    V«S«    W     <M=1> 
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RATIO    OF    CUT-OFF    SLOPES    OF    MTBC    AND    MB 
FUNCTIONS    V»S«    H    <M=1 > 
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Figure    23    -      RATIO    OF    CUT-OFF    SLOPES    OF    MTBC    AND    MB 
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RATIO    OF    THE    C'JT-OFF    SLOPES    OF    MTBC    AND    C 
FUNCTIONS    V«S«    W    <M=1> 
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Figure    25    -      RATIO    OF    THE    CUT-OFF    SLOPES    OF    MTBC    AND    C 

FUNCTIONS    V«S«    W    <M=2> 
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Figure    26    -       THE    DIFFERENCE    BETWEEN    THE    STOP-BAND 
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Figure    27    -       THE    DIFFERENCE    BETWEEN    THE    STOP-BAND 
ATTENUATIONS    OF    MTBC    AND    C    FUNCTIONS    V»S»    tf     <M=2> 
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Figure  28  -   RATIO  OF  THE  CUT-OFF  SLOPES  OF  THE  MTBC  AND  MC 
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RATIO    OF    THE    COT-OFF    SLOPES    OF    THE    MTBC    AND    MC 
FUNCTIONS    V«S«    W    <M=2> 


59 


K  -C!i~.p! 


inF-n  i 


s  'N T 


i.  :  1  •-  V  :  I 

INC1-! 


f.O 


*.1 


1-H 


u3o  Cr^d./%e.Cj 


f-r 


5-3 


Figure  30  -   RATIO  OF  THE  CUT-OFF  SLOPES  OF  THE  MTBC  AND  TBC 

FUNCTIONS  7«S«  W  <M=1> 
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HATIO    0?    THE    C0T-0FF    SLOPES    OF    THE    MTBC    AND    TBC 
FUNCTIONS    V»S«    W    <tf=2> 
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Figure    32    -       THE    DIFFERENCE    BETWEEN    THE    STOP-BAND 
ATTENUATIONS    OF    MTBC    AND    TBC    FUNCTIONS    V»S«     H    <M  =  1> 
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Figure    33    -       THE    DIFFERENCE    BETWEEN    THE    STOP-BAND 
ATTENUATIONS    OF    MTBC    AND    TBC    FUNCTIONS    V«S«    H    <M=2> 
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IV.   COMPUTER  PROGRAM 


A.   INTRODUCTION 


in  the  previous  chapters,  expressions  for  MTBC  filter 
are  derived  and  the  performance  of  this  filter  is  compared 
with  various  filters. 

In  this  chapter,  A  computer  program  is  developed  to 
implement  MTEC  filter.  It  is  pointed  out  that  all  five 
filters,  which  are  used  to  compare  with  the  MTBC  filter,  are 
the  special  cases  of  the  MTBC  filter.  Thus  the  program 
developed  in  this  chapter  for  the  implementation  of  the  MTBC 
filter  may  also  be  used  to  implement  any  one  of  these  five 
filters. 

A  sample  problem  is  worked  out  to  illustrate  the  use  of 
the  plots  presented  in  earlier  chapters,  the  use  of  the 
computer  prcgram,  and  to  point  out  the  flexibility  offered 
by  the  MTBC  filter  to  the  filter  designer. 


B.   COMPUTER  PROGRAM 


In  the  most  general  form  the  transfer  function  of  MT3C 
filter  is  given  by  eg.  (2-27)  ,  which  is  repeated  here  for 
convenience . 

1=1 


F(ju>) 


f  Cw.1- 10*) z  •*  ICwM)  ulW  Cn\vO) 

1  =  1  i=i 


H 
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where  : 


4lt_ 
w  =  Location  of  l—  inserted  zero 


m  =  Order  of  inserted  zeros 

n  =  Order  of  filter 

k  =  Weighting  factor  of  Transitional  Butterworth 
Chebyshev  filters 

When  k-n0  and  m=0,  |F(jw)|2  is  identical  to  the 
Butterworth  function.  With  k=n  and  m*0,  |F(jw)l2  is 
identical  to  the  modified  Butterworth  function  which  is 
discussed  by  BUDAK  and  ROY  [1]  -  [2].  When  k=0  and  m=0, 
|F(jw) |2  is  identical  to  the  Chebyshev  function.  With  k=0 
and  m*0  Modified  Chebyshev  function  may  be  implemented  which 
is  discussed  by  AGARVAL  and  SEDRA  [3]. 

When  m=0  and  k#0,  |F(jw) |2  is  identical  to  the 
Transitional  Eutterworth-Chebyshev  filter.  Finally  allowing 
both  k  and  m  to  vary,  Modified  Transitional 
Butter worxh-Chebyshev  filters  may  be  implemented. 

The  program  consists  of  two  main  parts.  In  the  first 
part  the  analog  filter,  which  is  given  by  eg.  (2-27),  is 
implemented  and  the  frequency  response  is  plotted. 

In  the  second  part,  the  analog  transfer   function   F(jw) 
is   first   predistorted   then   transformed   into  z-domain  by 
algebraic    substitution     method      (using     Bilinear 
z-transf ormation)  to  obtain  H  (z)  .   Digital  transfer  function 
[  H  (z)  ]  is  then  factored  into  second  order  cascaded   stages. 
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Finally,  frequency  response  curves  (  |F(jw)|  v.s.  w,  20 
log  |F(jw)|  v.s.  v.  and   F  ( jw)  v.s.  we  )  are  drawn. 

3ilinear  Transformation  is  preferred  over  the  other 
available  algebraic  substitution  metnods  {i.e.  Impulse 
invariant,  Marched  z-transf ormation)  in  obtaining  H  (z)  , 
mainly  for  the  following  reasons  [4]  : 

(1)  it  has  the  property  that  realizable  stable 
continuous  systems  are  mapped  to  realizable  stable  digital 
filters. 

(2)  iiidehand  sharp  cut-off  continuous  filters  can  be 
mapped  to  wideband  sharp  cut-off  digital  filters  without  the 
aliasing  in  the  frequency  response. 

(3)  After  Bilinear  Transformation,  the  relation  between 
the  analog  and  digital  frequencies  is  given  by 

3y  choosing  cut-off  frequencies  approaching  to  f»/2  where  f& 
stands  for  sampling  frequency,  extremely  sharp  cut-off 
slopes  may  be  obtained.  As  an  example  plot  of  cut-off  slope 
of  bilinearly  transformed  Butterworth  filter  v.s.  cut-off 
frequency  is  given  in  PIG.  34. 


C.   3ZQUIRED  DATA  CAHDS 


The  data  cards  required  to  use   the   program   are   given 
below. 


Card  1  :  Values  of  n,m,k  in  312  format 
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CUT-OFF    SLOPE    OF    BILINEARLY    TRANSFORMED 
BUTTERWOETH    FILTER 
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Card  2  :  Location  of  inserted  zeros,   w;   where   i=1,2,...,m 
in  8F10. 5  format 

Card  3  :  Initial  and  final  values  of  the  frequency  to  be 
used  for  the  frequency  response  plot  in  2F10.5  format 

Card  4  :  Number  of  solutions  required  in  13  format 


D.   REQUIRED  SUBROUTINES/FUNCTIONS 


In   addition   to  the  built-in  subroutines,  the  following 
IBM  source  library  subroutines  are  used. 

1 .  POLRT 

2.  PLOTP 

3.  PSUB 

4.  PMPY 

5.  PADDM 


E.   DESIGN  EXAMPLE 


Suppose  we  want  to  design  a  digital  filter  with  a  cut 
off  slope  of  15  and  minumum  attenuation  in  the  stop  band  of 
60  Db. 

Possible  solutions  for  various  types  of  filters  are 
given  in  table  II. 
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TABLE  !!•  POSSIBLE  SOLUTIONS  OF  DESIGN  EXAMPLE 


Or 

n. 

m 

at. 

. 

£T0P-&#klD 
ATTENUATION 

B 

43 

— 

— 

/SI 

Ci reader  -than 

M6 

c*jO 

2 

Hi 

/o~.0 

\8o±& 

3 

-M4 

^tf.05" 

|20  dS 

c 

7 

/W5" 

MC 

7 

1 

122 

/f.o 

^</£ 

2 

J.2Z 

2(3.5" 

Je^s 

13C 

2 

/7"£T 

0     63  J8 

H1SC 

7 

2 

/f.4 

6o</8 

2 

1 

1 

60  J£ 

i 

i 

2            /.// 

£// 

60  d'S 
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The  following  input  data  specifies  a  HTBC  filter  of 
order  with  111=2,  w0=1.36,  k=1,  sampling  period  of  1  sec, 
100  solution  points,  and  frequency  response  plot  from  0  to  3 
Hz. 


I 

2,3,4,5 

s 

1  1  S  1   9  |  10 

II   1  12  |  13  |  14  |  15  |  16  |  17  |  18  |  19  |  20 

21  |22|23|24|2S|26|27|2S|29|  30 

?,      .4 

I 

1       1       1 

1        1        1       1       1        1       I       1        1 

' 

l      I      I 

h.xZ-.t 

11.11      ill.  .3,6 

I      I      I 

\O\.\o 

.I....  i3i.iO 

I 

0<O<    i 

1  1   1 

1    1    I    1    1    1    1    1    1 

I    I    I 

1  1   1 

1    I    1    1    1   1    1    1    . 

Frequency  responses  of  three  of  the   possible   solutions 
of  the  design  example  are  given  in  Fig. 35. 
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Figure  35  -   MAGNITUDE  RESPONSES  OF  (a)  MB  FILTER  N=8,  M=1f 
(b)  MC  FILTER  N=7,  M=1,  (c)  MB  FILTER  N  =  7,  M=2 
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v-       Hill    DOMAIN    RESPONSE    OF    DIGITAL    FILTERS 


A.   INTRODUCTION 


In  previous  chapters  we  have  investigated  the  problem  of 
desining  digital  filters  in  the  frequency  domain,  i.e.  to 
meet  given  frequency  domain  specifications.  A  filter 
designer  should  always  consider  the  transient  response 
characteristics  of  its  filter.  There  are  many  applications, 
such  as  digital  MTI  filter,  for  which  one  is  interested  in 
the  transient  responses  of  filters  that  are  specified  in  the 
frequency-domain.  Time-domain  and  frequency^ domain 
characteristics  of  a  filter  will  work  against  each  other. 
Filters  close  to  the  ideal  frequency  characteristic  can  be 
designed.  Filters  whose  time  characteristic  is  close  to  the 
ideal  can  also  be  designed,  but  filters  close  to  both  cannot 
[12]. 

In  this  chapter  time-domain  response  of  digital  filters 
will  be  discussed  and  it  will  be  shown  that  the  location  of 
the  transfer  functions  poles  has  a  profound  effect  on  the 
transient  response  of  the  filter. 


B.   TRANSFER  FUNCTIONS'  POLES  AND  TRANSIENT  RESPONSE 


Given   a   transfer   function   of   a   digital   filter   in 
factored  form 
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where  z,  *s  and  pf  '  s  are  the  zeros  and  poles,  respectively, 
of  the  filter.  In  general,  the  input  signal's  Z-transform 
is  of  the  form 


X(4>  = !lfe) 


H 


where  g; 's  are  the  poles  of  the  input  function.  The 
response  of  the  system  may  be  obtained  from  the  transfer 
function  relationship 


Y(*)=  Hfr)  ■)<{*)-  K. 


fe.p) . . .  U-?n) 


Hit) 


C*-q,)---(l-lO 


(r-3) 

For  simplicity,  we  assume  that  all  of  the  poles  of  the  Y (z) 
are  distinct,  and  making  partial  fraction  expansion  of  Y  (z) , 
we  obtain 


iM  =  k, 


teifc 


_] —  + 


K-i 


+ 


C,  * 


+  ...  + 


C** 


The  response  of  the  system  may  be  decomposed  into  [5]  two 
parts  called  the  input  signal  mode  (ISM)  and  the  system  mode 
(SM)  as 


Y«W- 


a-p, 


+ 


+  ••■  + 


Combining  eg.  (5-4)  with  eg.  (5-5)  and  eg.  (5-6) 


Eg.   (5-7)   indicates  that  the  response  of  any  linear  system 
to  any  input  will   contain   modes   generated   by   the   input 
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signals  poles  and  the  system  transfer  function  poles.  The 
time  domain  response  of  the  filter  may  then  be  found  by 
taking   the    inverse    Z-transform   of   eg. (5-7) 

y(k)  =  k  sk  +  ^,SMa)  *  ^  c\l)  (s-t) 

where 

Equations  (5-9)  and  (5-10)  indicate  that  a  fast  responding 
system  is  one  in  which  all  of  the  system  transfer  function 
poles,  p; ,  are  sufficiently  smaller  than  unity  in  magnitude, 
in  order  that  the  system  mode  will  decay  to  zero  rapidly. 
On  the  other  hand,  a  slowly  responding  system  is  one  in 
which  the  system  mode  decays  to  zero  very  slowly  (i.e.  at 
least  one  of  rhe  p   is  close  to  unity  in  magnitude) . 

Those  poles  closest  to  the  unit  circle  will  be  called 
the  system's  dominant  poles  because  they  tend  to  dominate 
the  characteristic  of  the  resultant  transient  response. 

The  dominant  poles'  response  property  is,  as  stated  by 
Cadzow  [14]  "The  response  time  of  a  linear  discrete  system 
is  directly  dependent  on  the  locations  of  the  system 
transfer  function's  dominant  poles.  Depending  on  the 
particular  response- time  requirement  for  a  given 
application,  we  then  have  to  correspondingly  locate  the 
dominant.  poles  of  the  transfer  function.  A  fast  responding 
system  necessitates  dominant  poles  of  magnitude  much  less 
than  one  ". 

The  poles  of  the  MTBC  filter  for  the  orders  3-11  and  for 
various  values  of  m  and  w0  are  given  in  Table  III. 
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TABLE  III*  POLES  OF  HTBC  FILTER 


1 

n 

i 
rn 
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J.         i 
i 
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TABLE    III*    CONTINUED 
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Although  this  partial  fraction  expansion  method  helps  us 
to  understand  the  importance  of  the  poles  of  the  systems 
transfer  function  in  transient  response  analysis  of  digital 
filters,  evaluation  of  residues  of  corresponding  poles  of 
partial  fraction  expansion  is  not  a  trivial  problem. 

We  believe  that  the  so  called  ■  transfer  matrix  ' 
method,  which  we  are  about  to  discuss,  is  more  suitable  for 
digital  computer  simulation. 

Given  a  digital  filter  weighting  seguence  h (n)  and  the 
input  seguence  x(n),  the  response  of  the  digital  filter  may 
be  obtained  by  convolution  summation 


rt. 


H 


Eg.     (5-11)    may    be    written   in    matrix   form   as 

r KM         o  o   •  •  •     © i ^m 

kto)   h(o)     o  •  •  •  o 


gfc) 


3^  j 


WW      WGi-0     kU)  •  •  •  W«) 


Y<0 


H 


Or 

i}-  =  fi.  .  X  (s--ft) 

where  G  is  the  systems  transfer  matrix,  which  is  defined  as 


-ww     o        o  ...    o- 

KUj   WW    o  •  •  ■   o 


G  = 


.  kit)     wM    M*-0 


WO 


H 


To  find  the   systems   response   using   eg.   (5-13) ,   systems 
transfer   matrix   must  be  available.   To  obtain  the  weigting 
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sequence  (impulse  response)  of  the  system  to  form  the  system 
transfer  matrix  in  terms  of  the  transfer  function's 
coefficients,  the  difference  equation  of  the  system  is 
solved  for  the  impulse  input,  i.e., 


m 
ku)  =  O.o  -  J_  hi   MCr\-L)         ,  \)o     (sws) 

where  a-  and  b{  are  the  numerator  and  denominator 
coefficients,  respectively,  of  the  transfer  function  of  the 
digital  filter. 

A  computer  program (  FORTRAN  )  is  developed  to 
investigate  the  transient  response  of  the  digital  filters, 
using  equations  (5-13)  and  (5-15).  Program  listing  is  given 
in  Appendix  3. 

Step  responses  of  MTBC  filter  for  various  values  of  n, 
m,  and  w0  are  given  in  figures  36-39.  Figure  36  and  37 
indicate  that  increasing  the  values  of  m  and  w0  also 
increases  rhe  overshoot  and  settling  time,  but  doesnt  have 
any  significant  effect  on  the  rise  time.  The  rise  time  tends 
to  increase  with  increasing  order  of  the  filter. 
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C.   S0HH1BI 


Transient  responses  of  digital  filters  depend  en  the 
position  of  the  poles  of  its  transfer  function.  1  fast 
responding  filter  has  poles  of  magnitude  much  less  than  one. 
Modification  cf  ail  pole  filters  increase  settling  time, 
decrease  peak  overshoot,  and  doesnt  significantly  affect  the 
rise  time  of  the  filter's  step  response. 


30 
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Figure    36    -      STEP    RESPONSE    OF    MTBC    FILTEB(N=7,     W    =1«U6) 
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Figure  37  -   STEP  RESPONSE  OF  THBC  FILTER  (  N  =  5,  M= 1 ,  K=1) 
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Figure  38  -   STEP  RESPONSE  OF  MTBC  FILTER  (  M=1,  W  =1«36) 
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Figure   39    - 


STEP    RESPONSES    OF    B,    MB,    MTBC    FILTEfiS  (    N=5, 
H=1,     W    =1*36) 
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APPENDIX  A 


COMPUTER  PROGRAM  LISTING 


C 
C 

C 
C 


C 

c 
c 


c 

c 
c 
c 


c 

c 
c 


: 
c 
c 


CIMENSICN  RNUMR(25  )  ,  RNUM I ( 25 ) , RDNUMR ( 25 >  ,RCMMI(25)  . 
*C1(25)  tC  (25  )  ,R1(25) ,F (25) 
CIMENSICN  X(50)»Z(50),W(50),Y(50) 
CI  MENS  I  C  Is  FCCTR(50) TRCOTI(5C),COF(5C ) 
CIMENSICN  CM50)  ,CY(5C) ,CZ(50) 
CIMENSICN   RCOTDR(50 >,RGOTCI(50> 
CIMENSICN  CYH50) 

CIMENSICN  FX(100),PY(1CC),RC(1C0),RC1<100) 
CIMENSICN  FN ( 25  I ,RC(25) 
OIMENSICN  XT  (25)  ,YT(25),AX(25),AY<25) 
CIMENSICN  CUMMY( 100),CLMMX(100)tITB(20)tRTB(30) 
CIMEISSICN  CLMKAY(ICO) 
CCHPLE)  Fl.PiCtCl 

FiAC  THE  CFCEF  CF  FILTER, CRCER  CF  INSERTED  ZERCS  AND  THE 
WEIGHTING  FACTCR 

RSAC (5  ,5C1  )  N,M,K 
5C1  FORMAT<2I2) 

PEAC  LCCA7ICN  CF  ZERCS 

FEAC(5,5CC>  (W(I),I=1,M) 
5CC  FCRMAT  (EPIC. 5) 

RE/O  INITIAL  ANC  FINAL  VALUES  CF  THE  FREQUENCY 
(10  BE  LSEC  IN  FREQUENCY  RESFCNSE  CALCULAT ICNS  ) 

PEAC (5  ,5C2  I  W8EGIN,V<LAST 
5C2  FCPMAT  (2F1C5) 

PEAC  THE  NLNE  EP  CF    SOLUTION  POINTS 

PEAD(5,5CS)  NPCINT 
5C<  FCRMATU2) 


PEAC  SAMPLING  PERIOD  CF  THE  DIGITAL  FILTER 

7 


REA0(5  .7C  )  T 
FORMAT  (F1C.5) 


XPCINT=NPCINT 
WCELTA=('*LAST-WBEGIN  >/XPOINT 

CAlCLLATICN  CF  NUMERATOR  FCLYNCMAL 

ICIMY=1 

Y(  1)  =  1.C 

IF(M.EC  .C  )    GC  TO  401 

CC  2  I=lt> 

X(  l)=Vi  (I  )**4.0 

X {3)--2  ,C*V  (I  )**2.C 

X(2) =C.C 

X(5)=1.C 

X(4)=C«C 

CALL  FfrFY  (2,IDIMZ,Y,ICIMY,Xt5) 

CC  3  J=1,ICIMZ 

Y( J) =Z  (J) 
2  CCMIME 

ICIMY=1CINZ 
2  CCNTIME 

CCNVERSICN  CF  INDEPENCENT  VARIABLE  FROM   fc   TO   S 

CO  26  1  =  2  ,ICIMY,4 
Yd  >=-l.C*Y(I  ) 
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26  CONTINUE 
C 

C   FINC  THE  F-CCTS  OF  NUMERATOR  SQUARED  FUNCTION 
C 

PMJM*XCII'>-1 

CALL  PCLR7  ( Y  ,C0P  ,  MNU N  ,ROCTR  ,RCOTI  , I ER ) 
C 
C 

C   SELECT  THE  RIGHT  HALF  PLANE  RCOTS  OF  NUMERATOR  SQUARED 
C   FUNCTION 


C 


L  =  l 

CO  2  0  I=1,MNUM 

IFCRCCTMI  ).GT.O.O)  GO  TO  20 

RNUMR(L)=FCC7R(I ) 

RNUMKL  l  =  PCCTI(I  > 

L  =  L+1 
2C    CONTINUE 

MN=MNUN/2 

LK  =  L 

LZ=L-1 

IF(LZ.EC.NN  )    GO    TO    2  2 

LL=1.C 

CO    35     I=1,NNUM 

LX=L+1 

LV=I  +  1 

IF (RNUM  KLD.EQ.ROOTK  I  )  .OR  .RNUM  I  (  LL  >  .EG.  (-1  .0*ROGTI 
MI)))     C  C    T  C    3  5 

IF(LL.EC.  1)    GO    TO    26 

LJ=LL-2 

IF(RNLMI  (LJ).EC.RCCTKI  ) .OR .RNUM  I  < L J  )  .  EG  .( -  1  .C*R0G7I 
MI  )  )  )  GC  TC  25 

IF (LL  -  E  C  .2  )    GO  TO  26 

LI =LL-4 

IF(RNUMI(UI).EQ.RCCTI (I  )  .0  F  .RNUM  I  (  LI  ).EG  .(-1  .0*R0CTI 
MI  ))  )  GC  TC  25 
It    RNUMR( L)  =FCCTR(I ) 

RNUMKL  >  =  RCCTIU  ) 

RKUMBCLXJ =FCCTR(UY) 

RNUMI(L>)=FCCTI(LY) 

IF{LX.EdvN)    GO    TO    2  2 

LL=LL+2 
25    CONTINUE 
2  2    CC    2  2     JM,MN 

R€2  )=C>FL>  (FNUMRU  )  ,BMJNI  (I)  ) 
22    CONTINUE 

CAUL    MKFCL     (MNtRfC) 

^1  =  ^+1 

C(MN1)=CMFUX(1. 0,0.0) 

CC    31     1  =  1,  CM 

RCCI ) =P  E AL ( C( I ) ) 
21    CONTINUE 

GC    TC    5 
4C1    RC(1)=1.C 

MMM 

NN  =  C 

RNUMR(  1)=C.C 

RNUMH  1  )=0  ,0 

C^lCULATICN    CF    DENOMINATCR    POLYNOMIAL 

5    NDNUM=2*N 

IF(M  .EC  .C  .4NC.K.EC.N )    GO    TC    402 

CCNST=1.C 

IFCH.EC.C  )    GO    TO    7 

CC     11     1  =  1, M 

CM  L7MV»(I)**2. 0-1.0)  **2.0 

C0NS7=CML1*C0NS7 
11    CONTINUE 
7    IF(K.EC.C)    GC    70    52 

KK=2*K 

KKK  =  KM1 
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c 
c 
c 
c 


c 

c 
c 

c 
c 
c 


c 
c 

c 


c 
c 

: 
c 


c 

c 
c 


c 
: 
c 


CM(KKK)=CCNST 
CG  12  1  =  1,  KK 
CMI  )=C.C 
12  CONTINUE 

GO  TC  52 
52  CM( 1)=1.C 

KK=0.C 

KKK=1.C 
52  NK=N-K 

NKKl=NK4l 

DETAIN  THE  CFEEYSHEV  FCLYNCM I A L , SQUARE  IT  AND  NULTIPLY 
THE  RESULT  EY  EUTTERkQRTF  SQUARED  FUNCTION 

CALL  OeSV(CYlfNK) 

CALL  FNPY  (CY,NKK,CY1,NKK1,CY1,NKK1) 
CALL  FNF>  <CZ,IDIMCZ,CY,NKK,CM,KKK) 

CETAIN  CENCMINATOR  SGUAPEO  FUNCTION 

CALL  PiBCCf CZtIDIMZ,Y,IDIMY,1.0,CZ,IDIMCZ) 

CCNVERSICN  CF  INOEPENCENT  VARIABLE  FROM  W  TC  S 

t    CG  17  I*3,ICIMZ,4 
Z(I)=-1.C*Z(I) 
17  CCNTIME 
6  NONUMICIN2-1 

FINC  THE  FCCTS  CF  DENOMINATOR  SQUARED  FUNCTION 

CALL  PCLRT  ( Z, COF  ,  NDNUM , ROOTCR  ,ROOTD I  ,  I  ER ) 
GC  TC  4C2 
4C2  CALL  RCCT  (  NCNUM, RCCTCP  ,RGQTCI ) 


SELECT  THE  FIGHT 
FLNCTION 


HALF  FLANE  RCOTS  OF  DENOMINATOR  SQUARED 


AC2    J  =  l 

CO    21     I=1,NCNUM 

IF(RCCTCF  (  I  )  •GT.O.O)    GO    TO    21 

RDNUMF  (J  )=FCCTCR(I) 

RDNUMI  < J)=RCCTDI  <  I  ) 

J  =  Jf  1 

21  CCNTIME 
NC=NCMN/2 
CC    22     1  =  1, NC 

fild  )  =  CNFLXRDNUMR(I  )  ,RDNUM  (I  )  ) 

22  CONTINLE 

CALL    NXKFCL     (MC,R1,C1) 

NCl=MC4l 

CKMC1  l=CVFLX(l,C,C.C) 
47    FCPMAT  (2>  ,25F11.3//) 

CO    32    1  =  1, NCI 

RCK  I  )=PEAL(C1(I  )  ) 
32  CONTINLE 

NCRMALIZATICN  CF  NUMERATOR  POLYNOMIAL 

FACTOR=RCl(l)/RC( 1) 
CC  ICC  1=1,  CM 
RC(I  )=PC  (I  )*FACTGR 
IOC  CONTINLE 

WRI T  E  <  6  ,57C) 

C  l'f  PUT  SECT  ICN 

57C  FCRMAT  (  '  1'////) 

U  R  T  T  C  f  £    C  6  C  ) 

5tQ    FQRMATC2X,'     LOW-PASS    PROTOTYPE  (CONT I NUOUS )     FILTER'//) 

bRITEU  ,571) 
571     F0RM*T(2X,'    CRDER    OF    FILTER     •//) 
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c 
c 
c 


1 


6C 
57 

57 

57 
8C 
57 

57 

51 
51 

51 
52 

51 


521 
576 
51C 
577 


FO 

V»R 
FC 
WR 

-  : 

*  = 

WR 

FC 

FC 

-  : 
FC 
( 

UR 
WR 
FO 
( 

kR 
'*R 
FO 
*R 
FC 
CC 
WR 
FC 
CC 
kH 
FC 
'*R 
CC 
l»R 
CC 
bP 
FC 
WR 
FC 

*R 

c  - 


IT5( 
RMAT 
ITE  < 
RMAT 
ITE  ( 
IT5( 
RMAT 
ITE  ( 
ITE( 
RMAT 
ITE  { 
RMAT 
ITE  ( 
R^AT 

:esc 

ITE( 
ITE( 

RMiT 

:zsc 

ITE  ( 

ITE( 
RMAT 
ITE  ( 
R^AT 

52C 
ITE( 
RMAT 
NT  IN 
ITE  ( 
RMAT 
ITE( 

521 
ITE( 
NTIN 
ITE  ( 
R'<AT 
ITE< 
RMAT 
ITE  ( 
RMAT 
ITE  ( 


6,6 
(2X 
t  ,5 
(2X 
6  ,6 
6  ,5 
(2X 
6,6 
6,5 
(2> 
6  ,6 
(6F 
6,5 
<2X 
ENC 
6,4 
6  ,5 
(2X 
ENC 
6,4 
6,5 
<2X 
6,5 

(ex 
1  = 

6,5 
(5X 
LE 
6,5 
<2X 
6,5 
1  = 
6,5 
LE 
6,5 
12X 
6,5 
<2X 
6,5 
(2X 
6  ,5 


CC)  N 
,  12//) 
72) 

,  ■  CRQER  CF  INSERTED  ZEROS  •//) 

CC)  M 

72) 

,  »  WEIGHTING  FACTOR  CF  TBC  FILTER  '//) 

CC)  K 

74) 

t  ■  LOCATIONS  CF  ZEROS  '//  ) 

CI)  (w(I  )  ,1=1, M 

1C.5//) 

76  ) 

,'  COEFFICIENTS  CF  NUMERATOR  POLYNOMIAL 

ING  CROER  )     •//) 

7)  (RC( I  )  ,I  =  l,yNl) 

75) 

,  •    COEFFICIENTS    OF    DENOMINATOR    POLYNOMIAL 

If^G    CRDER     )     '//) 

7)     (RCKI  )  ,1  =  1, MCI) 

11) 

,'RCOTS  CF  NLVERATCR'// ) 

12) 

,'REAL  PART  •  ,4-X, • IMAGINARY  PART  '//) 

1  ,¥h 

12  )  RNUMRd  ),RNUVI(  I  ) 

,F10.5,4X  ,F1C5/  ) 

14  ) 

,'RCCTS  CF  CENCMINATCR'// ) 

12  ) 

It  PC 

12)  RDNUMR(I) tRCNLMI(I) 

76) 

,  •  INITIAL  ^ALLE  CF  FREQUENCY  '//) 

1C)  WBEGIN 

,F1C.5//1 

77) 

,«  FINAL  VALUE  CF  FREQUENCY  •//) 
1C)  «LAST 


CUTTING  SECTICN 

WRITE(6  ,5<C) 
55  C  FORMAT  ('l1) 

CALL  PPLCT  TC  FIND  FRECLENCY  RESPONSE  VALUES  CF  THE 
FILTER 

CALL  PFLCT(RC,RCl,W5EGIN,wCELTA,NPOINT,PX,PY,fDl,MNl) 
CALL  FLCTF  ( FX , PY ,NFC  INT , C  ) 

PScCISTCRTICN  FOR  BILINEAR  TRANSFORMATION 

SCALE=7A^  (T/2. ) 
CO  777  LK=1,MN1 
EXF=LK-1 

RN(LO=RC(LK)/(SCALE**EXP) 
777  CCNTINLE 

CO  776  LK=1,MD1 

EXP=LK-1 

R0(LK)=BC1(LK)/(SCALE**EXP] 

77c  CCNTINLE 

n  t  =  m  fv 

CALL  ZCMN  (RN  ,RD  ,NK  ,f«T,WBEGIN,WLAST,  kCELTA  ,NFCINT,7) 

CALL  XFCRMRC,RC1,1,XT,YT,*N,MC,1CCC.C,:2.C) 

CALL  PFLCT(XT,YT,5C0-C,10.C,1CC,AX,AY,NC1,VM) 

CALL  FLCTF  (AX, AY, ICO, C) 
STOP 
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c 

c 

: 
: 
: 
c 


c 
c 

r 


C> 


ENC 


SL3R0LTINE  OBSV  <CY,NC) 

5t?59JrIlt?  JC  FIND  Tl~6  CCEFFICIEM5  OF  CFEEYSFEV 

FClYNC*IAL  CF  GIVEN  CPCER 

NC  :  CRCEF  CF  CFEEYSFEV  FCLYNOIAL 

:  CALCILATEC  COEFFICIENTS  (CESCENOING  C^CER) 

CIVENSICN  CX(50)  ,CY(5C),YY(50),ZZ(50  ),Z(50) 

N  N  =  N  C  4  1 

CX(1)  =  1  ,c 

IF(NC.EC.C)  GG  TC  2 

CY(1)=C.C 

CY(2  1*1. C 

Ic ( NC  .  EC  .  1  )  GG  TC  19 

YY (1  )  =  C  .C 

YY(2)=2.C 

CO  5  I  =  2  t NC 

11=1+1 

11=1-1 

CALL  ?NFY(Z,IDIMZ,YY  ,2.CY,I  ) 

I C  I M  Z  Z  =  I  C  I  fi  Z 


c  o 


1  •  ICINZ 

) 


ZZ(J)=Z( 
CONTINLE 

CALL  FSia(Z,I3IMZ,ZZ,i:iVZZ.CX,U) 


6 
LFLATE  FCLYNCMALS 


CO  7  J  =  1  , 1 
CX(J  )  =  CY ( j  I 

7  CONTINLE 

CG  8  J=l  ,  ICI vz 
CYU>  =  Z<JI 

6  CCNTIME 

5  CONTINLE 
l<  RETURN 

i  CY(1)=C> (1 ) 
GC  TC  IS 
ENC 


SLi 
FIL 

PCI 
PC 
N  : 
V  : 
Fl 
F2 
T  : 


SLSPCL7INE 
TC 


ZDMN(RCtRCl,N,M,FltF2,WDELTA,NP ,x ) 

FIND  ANC  FLCT  FREQUENCY  -RESPONSE  CF  DIGITAL 


ROUTINE 
TEPS 

:  COEFFICIENTS  OF  DENCPINATGR  OF  CONTINUOUS  FUNCTION 
COEFFIC  I6NTS  OF  NUVERATGR  CF  CONTINUOUS  FUNCTION 

CPCEP  CF  CEN0WIN4TCR 

CPCER  CF  NUMERATOR 
:  INITIAL  FREQUENCY 
:  =  INAL  FPECLENCY 

SAVFLINC-  PERIOD 


CI 
CI 

PC 
X° 

ex 

CI 

CI 
FF 

:i 
,p 
ci 

:: 

CO 


VENSICN  ECC(25) 

MENS  ION  FC(25),PCH25>tZX(25),ZXl(25)  ,  XJ(25),YJ(25>, 

CTPC (25)  tPCCTPN  (25  )  ,RCCTIC(25) ,RGOTIN ( 25 ) ,  CY(25) , 

AY(25 ) ,>FA>DA( 25) , 

(25  »  ,YFJY(2  )  ,YP£YDM2  ) 

KENSICN  ACC(IO)  ,Aii(10),*22(  10 ),B11 (iC)  , 522(10) 

MENSICN  >AXIS(l5C0),YAXISH150C)tYAXIS2(l5C0), 

(15CC  )  ,CEL(  150C) 

MENSICN  >CCF  (3)  ,  X2CCF(3)  ,  COF  (  3  ) ,  RCCT  .<  (  2  )  ,  FCCT  I  (  2  ) 

C0T1  (2  )  ,;CCT212 ) 

VENSICN  CLfKAY(lCC) 

UBLE  PRECISION  W,Y 

M*GN/FAW/W(100) ,Y( 10C>»ICALL 

AL  K/Gh  ,  MGNH,NL.*i  ,NUN2 


39 


CX  ,C>,XJ , VJ 

NLN,CENQV; 


9C 


COMPLEX 

CCNPLEX 

MJ=M+1 

NJ=N+1 

NNN=N-> 

call  ez>ffm,rc,m,nvm  ,zx,1) 
call   s2xfpmrc1,n,nmn, zx1,c> 

sfact  =  zx<n,j)/zxi(nj) 

CC  9C  1  =  1, NJ 

ZXH  I  )=ZX1(I)*SFACT 

CCNTINLE 


WRITEU  ,6C) 
C 

C  COEFFICIENTS  CF  DENOMINATOR  PCLYNOPIAL 
C 


i: 

6C 

7C 


61 

62 

62 
64 


WRIT 
WRIT 

CALL 
CALL 
CALL 
CALL 
CALL 
F  =  F1 
CO  1 
X  =  F* 
VAGh 
FF3  = 
CG  6 
AC  =  A 
A1  =  A 
A2=A 
S  1  =  3 

E^  =  e 

NLMi 
NUM2 
NLV  = 
CeNO 
CENO 
CENC 
VAGh 
FH1  = 
PH2  = 
F!-2  = 
CC.NT 
VAGN 
PF(I 
FK  =  P 
XAXI 
YAXI 
>AXI 
MI) 
Y(  I  ) 
F  =  F  + 
CCNT 
WRIT 
FORM 
FGRM 
hRIT 
FORM 
(PRC 
WP-IT 
FORM 
WFIT 
*RIT 
FORM 
WRIT 
WRIT 
FOR* 
WRIT 
FGRM 


£(6  ,45  )  (ZXH  I  »  ,  I  =  1  ,NJ  ) 
EU  t6C) 

PCLRT(ZX,CCF,N, FCOTRN  ,RGOTIN,IER) 
FCL=T(ZX1,CCF,N ,RCCTRC,RCCTIO, I  ER  ) 
CCE^EN  (N,MX) 

F^CTCR (ROGTRN,RCCTIN,N,VX,AGO,A11  ,A22  > 
F4CTCR ( ROGTR C, P CCTI C, N,MX, 3 22, Bli, ECO) 

C  1  =  1 ,  NF 

T 

*1.C 

c.c 

c  c  T 

n  ( 

22  ( 
11  ( 
22  ( 
=  (A 
=  U 
CVP 
Ml  = 
N2  = 
N=C 
=  AE 
(AT 
(AT 
FF2 
INL 
F  =  2 
)  =  F 
h(l 
S(  I 
Sl( 
S2( 
=  XA 
=  W 
WCE 
INL 
EU 
ATI 
AT  ( 
E  U 
AT  ( 
TCT 
E  (6 
AT( 
E  (6 
EU 
AT  ( 
E(  t 
EU 
AT  ( 
EU 
AT  ( 


1,N 

J  ) 

J  ) 

J) 

J  ) 

J) 

2  +  1 

C-A 

LX  ( 

(E2 

(  1, 

NFL 

S(N 

if\  ( 

AN  ( 

+  (  F 

E 

C.« 

F2 

) 

)=F 

I  )  = 

I  )  = 

X  IS 

XIS 

LT  A 

E 

,6C 

■  1  ' 

2X, 

,7C 

2> 

YF 

,6 

2X 

,4 

t  62 

2X, 

,45 

,62 

2X  , 

U« 

IX, 


n*c:s(  xj+ai 

2)*SIN(XI 

NLM ,NUN2  ) 

♦i.0)*ccs(x)+ei 

-E2  >*SIN(X  ) 

X(CENGV1,CENCV2) 
AGh*CA3S(NL^/DENGV)  ) 
NLN2/NUM1))*57.2957  7951 
CENCX2/CENCM)  )* 57  .29  5  7795  1 
F1-PH2  ) 

ALCGIO(NAGF) 


NAGh 
f-AGNH 
(I  ) 
1(1) 


) 
) 

1CF12.3/// > 
) 


3ILINEARLY  TRANSFORMED  LCW-PASS 
DIGITAL  FILTER'//) 


COEFFICIENTS  OF  NUMERATOR 
(ZX(I  )  ,I  =  i,NJ) 


PCLVNCWIAL  »//) 


COEFFICIENTS  OF  DENOMINATOR  POLYNOMIAL  '//) 
(ZXK  I  )  ,1=1, NJ) 

NUMBER  CF  CASCADED  STAGES'//) 
NX 
2) 
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r 

c 
c 
c 


WRITEU  ,61) 
81  PgRMAT  <2X,  'COEFFICIENTS  OF  CASCADED  STAGES1//) 

.-    ^|l|^2e2)JG0(I),All(I),A22(I),Bll(I)te22<I) 

o:     rG«Mn  (i)<  t;pl2.5//  ) 

62  CCMIME 
WRITEU  ,63  ) 

63  FORMAT  (2X  ,  'ROOTS     CF    NUMERATOR '// ) 
CC    66    1=1 ,N 

WRITE<6,43)    POGTRNd  ),R0GTIN(I  J 

66  CONTINUE 
WRITEU  ,67) 

67  FORMAT <2Xf'  ROOTS  OF  DENOMINATOR  '//) 

WRIT  E  (  6  t  5  5  )  ROOTROd  )  ,RCCTIC(I  ) 

68  CCNTINLE 
WRITE (6  ,S£  ) 

96  FCRMAT(2>,'  SAMPLING  FRECUENCY  '//) 

WRITE(6,cc  )T 
9<E    FORMAT  (2X  ,F1C. 5//) 

N=NP 

WRITE(6,12  I 
12    FORMAT  (  •  1*  J 

CALu     FLCTF OAXIS  ,YAXIS1,\,C) 

WRITE(6,31  ) 

31  F0RMAT(//4CX,  •  ABS.  GAIN  V.S.  FRECUENCY') 
*RITE< 6,16) 

16  FORMAT  (  i  i«  ) 

CALL  FLCTF OAXIS,  YAXIS2,N,C) 
"RITE  (6  ,3  2  ) 

32  FCRMAT<//4CX,»  GAIN  (C3  )  V.S.  FREQUENCY  •) 
*RITE(6  ,14) 

14  FORMAT  (  '  1  •  ) 

CALL  PLCTP  (XAXIS  ,FI-,N,C) 
WRITE( c,32  ) 
52  FORMAT (//4CX,»  PHASE  V.S.  FRECUENCY  •) 
RETLRN 
ENC 


SUERCUTINE  CCEVEN(N,M ) 

SUBROUTINE  TO  FIND  THE  NUMBER  CF  CASCACEC  S' 
N:  Zmtz    CF  NUMERATOR  FCLYNGMIAL 
M:  CEGREE  CF  CENQMINATCP  FCLYNCMIAL 


'AGES 


NN=(N/2)*2 
I F  ( N  .  N  c  .  V  N  ) 
M=N/2 

GO  TC  2 
*=<N+1  )/2 

RETURN 
ENC 


10  TC  1 


SU 

CI 

*  ,Y 

3  ,X 

SL3R3 

FClYN 

*  :  FC 

»  :  CR 

N  f  M  =  N 

AND  C 

ZX:  R 

I=mCT 

i  fact 


BRCL 

MENS 

:jm 
fact 

UTIN 
CM  I  A 
LYNC 

CER 

-M   ( 

ENCM 
ESUL 
=  1 
=  C 


TINE 

ION 
)  ,2 
23  ) 


EZXFRM (X, N  ,NMM,ZX,  I  FACT  ) 
>(25),Y(25),ZX(23),XCUM(2) 
1125),Z2(23),X1( 25>,X2(25> 
,ZXi(23) 

FIND  ThE  BILINEAR  Z-TRANS FORM 


OF  A  GIVEN 


IAL 

CF  F 

TM 

IN  A  7 
TANT 

FCR 
FOR 


CLYNOHIAL 

:ic=ERENCE  EET^EEN 
CR 

FCLYNOMIAL 
NLVERATCP 
CE NOMINA  TCP 


ThE  CRCERS    OF  NUMERATOR 


M=V  +  NM+1 
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N2=M*1 

CO    1    1=1  .M> 

ZX(I  )=C.C 

CONTIME 

XDLN(1)=-1.C 


XDUMU 
YDIM1 
YOLM(2 
IZX  =  2 

CO  2 
CALL 
CALL 
CALL 


)  =  1 
)  =  1 
)  =  1 


,11) 

12) 

tI2) 


11  =  1  ,N2 

PCLE>F(XDUM,2,Hf,Xl 

FCLEXF(YCUM,2,N  ,X2, 

FNPXZltlZi  XI,  II, X2 
FACT  =  X  (II  ) 
CALL    FADCN  ( 22 ,1 Z2 , ZX , I ZX , F JCT  ,  Zl  ,  I  Z  ) 
CC    4    K=  1  ,  1  22 
ZX(K)  =  22(K  ) 
C3NTINIE 
NM  =  MM1 
N=P-1 
IZX=IZ2 
CCNTIME 

IF(IFACT.EC.C)     GC 
CALL    FCLEXF(YDUM, 
CALL     FNFY<2X1,MZ, 
CC    6    I  =  1 , N  2 
ZX(  I  )=ZX1<  I  ) 
CONTIME 
V=M2-1 
RETURN 
ENC 


TC    5 
2,NPM,XFACT,NXFACT) 
XFACT,NXFACT,ZX,IZ2) 


,PIN 
(25) 


SLERCL7INE     FACTOR(RRN 
CIMENSICN    RRN( 25)  ,RIN 
*    C(25),C1(25),AOO(25),A11(25),A22(25) 
CCNPLEX     C,C1 

IMTIALIZE    CCLNTERS 

IM  =  C 

IY  =  1 

K  =  0 

M=0 

IX  =  1 

NXR  =  C 

NXI  =  C 

IFtV.NE.C  )    C-C    TO     18 

CC    17    1  =  1  ,>> 

AJO<  I  )  =  1.C 

All( I)=C  .C 

A22(I  )=C.C 
17    CCNTINLE 

GC    TC    2C 
IE    IF(M.NE.l) 

AOO(  1>  =  1  .C 

A1K  1)  =-1 

A22(  1)=C. 

K  =  2 

NCIF=M> 

NXK=1 

GC    TC    1C 
U    CO    1    1  =  1  , 

IX1  =  IX-»1 

IFdX.GT. 

IF(RRN(IX).EC.RRN(IX1).AND.RIN(IX).EC 
* (-l.C*F  IN  (  1X1)  )  ) 

IMI»+1 

XR(IN)=-1.C*RRN( IX) 

IX=IX+1 


K,MX,A22,A11,A00) 
,  XR<25)  ,XI(25),YR(25),YI(25), 


GC    TO    16 
C*RRN< 1) 


ICC 

M    GC    TC    7 

)  .EC.RRN(IXl)  .AND.RIN(IX) 

GC    TC    2 
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11 

1C 


NXR=NX 
GC    TC 

iyi=iy 

YR(IY) 

YKIY) 

>R(IY1 

YKIY1 

IX=IX1 

M=NI  + 

NXI=NX 

IY=I Y 1 

CONTIN 

NT=NXF 

IF(NT. 

IF(MX. 

NCI(=  =  N 

IC  =  1 

CO    4-    I 

IC1=IC 

ACG( I  ) 

411  (  I  ) 

A22(  I) 

IC=IC+ 

CCNTIN 

N0IF1= 

NCCF=N 

NLAST= 

CC    5    I 

NCGF=N 

ACC(NC 

A1KNC 

A22(NC 

CCNTIN 

N6EGIN 

NCOMF= 

CC    6    I 

NC0MF1 

C(l)-C 

C(2) =C 

CALL    V 

ACCtNe 

All( NB 

422(NE 

NCGMF= 

CCMIN 

GC    TC 

NCIF=f 

IF(NXR 

CC    9    I 

ACCl  I  ) 

Alien 

422(  I  ) 

K  =  K+1 

CCNTIN 

IF(NXI 

NXK=NX 

CC     11 

NXK=NX 

IK=I+1 

C(1)=C 

D(2)*C 

CAuL    M 

ACC(NX 

A1K  NX 

A22(NX 

K  =  K  +  1 

CCNTIN 

CG    12 

NXK=NX 

ACCtNX 

A1HNX 

A22(NX 


P  +  l 

1 

+  1 

=RRN(IX> 

=  FIN  (IX) 

)=RRN(IX1) 

)=R1N(IX1) 

jl 

2 

1+1 

+  1 

LE 

■»NXI 

EC. NX)    GG 

CT.NT  )    GO 

T-NX 


NCIF 


TC 
TO 


F (IC  l+XR ( IC1  ) 
R (ICJ*XR(XC1) 


IF-f  1 

F*2 

-NXI 

CIF1,NLAST 

F4  1 

i  )  =  1.0 

1 )=XR(NCCF) 
1  )  =  C.O 

LAST+1 

BEGIN  ,MX 

CCNF+1 

LX(YR(NCOMP  )  ,YI(NCCMP) ) 

LX(YR(NCOMPl) ,YI(NCCMP1) ) 

PCL(2tD,Dll 

in  )=i.o 

IN)  =  FEAL(C1(2)  ) 
IN  )=REAL(C1(  1 
CNF1+1 


I  ) 


=  1 

41 
=  1 

=x 
=  x 
1 

LE 

NC 

CI 

NX 

-N 

CC 

IF 

IF 

IF 

LE 

=  N 

1 

=  N 

=  N 

NF 

NF 

*K 

EG 

EG 

EG 

NC 

LE 

2C 

X- 

.E 

=  1 

=  1 

=  X 

=  C 

LE 

.EC  .C)    GO    TC    IC 

r 

1  =  1, NXI 
K  +  l 

NPLX  (YR(I),VI(II) 

NPLX  (YF(IK)  ,YI  (IK) ) 

JKPCL     (2,0, CI) 

K  >  =  1  .C 

K)  =FEAL(D1 (2) ) 

K )=FEAL(01(  1)  ) 

LE 
I=K,NCIF 

i<  +  ; 

K)=1.C 
K  l=C  .C 
K  )=C.C 


NT 
CO 

,NXP 
.C 

R  (  I) 
.C 


GO    TC    15 
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12  CCNTIME 
QC    TC    2C 

2    IX  =  1 

LZ    13     1  =  1  ,NXP 
ACG( 11*1. C 
All! I  )  =  XR  CI  ) 
422(1  J=C.C 

13  CCNTIME 
NXX=NXF-H 

CC    14    I=NXX,vx 

I>l=IX4l 

C (1) =C*f LX(YR( IX  )  ,YI  (IX)  ) 
C(2)  =C*PL>  (YFCIX1  )  ,YI  (1X1)  ) 

CALL    M4KFCL     (2,0, CD 
4CC(NX>  1=1 .C 
A1K  NX>  >  =FEAL(C1  (2)  ) 
A22(NX>  )=FEAKCi(  1)  ) 
IX  =  IXH1 

njx=nx>+i 
14  continue 

2C  RETURN 

END 


C 

c 
c 

r 
C 

: 

L 

c 


SUBROUTINE  PQLEXP(  XX  ,IOIFXX,H, YY  ,1 DI  I»YY1 

SLSR.CLTINE  TC  FIND  THE  CC*=RS  CF  GIVEN  POLYNOMIAL 

XX  :  PCLYhCPIAL 

ICi  VXX  :  CIMNSICN  CF  XX 

V:  POKES  TC  EE  SAISEC 

YY: RESULTANT  POLYNOMIAL 

ICiNYY:  CIPEhSIQN  CF  RESULTANT  PCLYNOMIAL 

CIVENS1CN  XX(25)  ,YY(25  )  ,ZZ(25) 

IF  IN  .EC.  11  C-C  T0  4 

ICIMYY-1 

YY (1  )  =  1  .C 

IF (M. EC. CI  C-C  TC  2 

CC  1  1*1 fP 

CALL  FNFY  (ZZ,ICIfZZ,YY,  ICINYY, XX, IDT-'XX) 
CO  2  J  =  1  , 1  C  INZZ 
YY( J  I  =  Z  2  IJ  ) 

2  CCNTIME 
ICI^YY^ICINZZ 

1    CCMIME 

GC    TC    3 
A    CO    5    1  =  1  ,ICI^XX 

YY ( I  )  =  XX  (  I  ) 

i  ccntime 
ici.myy  =  i::*xx 

3  RETURN 
END 


SUBROUTINE    ccLOT(RC,RCl,weEGIN,«CELTA,NPCINT,?X,PY, 
*VC1,NN  1  » 
CI^ENSICN    FC(IOO)  tRCKlCC) 
CINENSICN     FX(lJG),PY(i:C) 
COPLEX    XMMtCNUM 
fc=*6EGIN 


CC    666     I=1,NPCINT 
b»h+fcCELTJ 

INITIALIZATION    CF    REAL    ANC    COMPLEX    PARTS    Jr    NUMERATOR 
A  NO    DENCHNATCR 

XNUMA=C.C 

x  n  l  y  e  =  c .  c 

XNUMC=C.C 
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_ 

>NLMD=C.C 

c 

XCNUNA=C.C 

XDMJNB=C.C 

XDNUMOC.C 

XCNUNC=C.C 

c 

c 
c 

INITIALIZATION    GF  CCLMERS 

NNLMC=1 

PNUM1=2 

NNUM2=2 

MNLM3=« 

c 

*»cnunc  =  i 

NCNUM1=2 

NCNUP2=2 

NDNUN3=4 

c 

XC=MNLNC-1 

X1  =  MNLM-1 

X2  =  MMN2-1 

X3=M NLN2-1 

XC0  =  NCMNC-1 
XC1^CMM-1 
X02=*DM*2-1 

XC3=*CNLN2-1 

OG    665    J=1,1C0 

IF(,VNLf-3.LE.NNl) 

IF(MMJf2.LE.MNl) 
IF(MNUM  .LE  .MN1) 
IF(MNLNC  .LE.NN1) 


I F  (  I*  C  N  L  I*  3 
IFCMDMN2 
IF(VCNLM1 
IFiVCMVC 


LE.MD1  ) 
LE.^Oi) 
LE.M01) 
LE.MD1) 


XNUNC=XNUNC+RC(MNUM3)*W**X3 
>NLNC  =  XNLNC  +  RC(MNUN»2)*W**X2 
XNUNE=XNUV8+RC(MNUM  1I*W**X1 
XNUNA=XMJNA+RC(MNUMC)*W**XC 

XDNUMC=XCNUMD+RC1(M0NLM3  )*fc**XD3 
XCMNC  =  XCMJMC  +  RC1(PCNUM2  )*W**X02 
XCMjNB=XCNUMe  +  RCl(^CMJNl  )*U**XD1 
XONUMA-XCNUMA+RCK  NCNIVO  )*U**XD0 


INCREMENT    CCLMERS 

nncmc^mk+4 

NMJMl  =  NNLNl+4 

MNUM2=^NL*2+4 

fCNUNG  =  »-CMN0+4 
yDNUMl  =  MCMM  +  A 
VCMIN2  =  NCMV2+4 
►0NU»'3»»'CNLI»3+4 

XC  =  XC*4  ,C 
X1=X1+4.C 
X2=X2*i .C 
X3=X3+4.C 

XC0-XCC*4,C 

>C1=XC1+4.C 
>C2=XC2+4.C 
XC2=XC2+4 .C 


C 
C 
C 


66 


IF(  (  NDM^C.GT.MDU  .ANC.  (MNLMO  .GT  .MNi  )  )     GG    TC    668 
CGNTIME 


CAuCLLATE    NAGMTUDE    CF    THE    TRANSFER    FUNCTION 

666  XMJM10MN4-XNUMC 
>NLM2=>NLNE-XNUM0 
XNUM=CFFLX(XNUM1,XNUM2I 

CNLN1  =  >DMNA-XDNUNC 
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ttt 


CNUM2=>CNLME-XDNUMD 

CNUM  =  ^FLX  (CNUM1  ,CNUN2) 

XMAGH=CAES(XNUM/DNLM) 

PX(I  \  =  V 

PY(I  )=2C.C*AL0G10(XMAGH) 

CONTINLE 

RETURN 

END 


C 
C 

C 

c 

c 
c 

r 

c 

c 
c 
c 
c 
c 
c 


SIB 

MY 
N7Y 
N1Y 
MY 
N1Y 
X 

\ 

w 

N 

XT, 


SL3RGL7INE  XFGRM  (  X  ,  Y  ,  MYPE  ,XT  ,  YT  ,  M  ,  N  ,  Fl ,  F2  ) 
RGUTINE  FCF  FREQUENCY  TRANSFORMATIONS 


c 

11 


2C 


FE  =  3 

PE  =  i 

PE  =  2 

PE  =  3 

PE  =  4 

COE 

CCE 

CRD 

GRD 

YT  : 

CIME 

NN=N 

N1=N 

M  =  M 

IF(N 

CO  1 

XM=N 

>T(! 

CCNT 

CC  2 

>N=N 

YT(I 

CCNT 

GO  T 

IF(N 

XF=M 

YF=N 

FACT 

K  =  M1 

rO  3 

XN=M 

K  =  K- 

XT(I 

CCNT 

K  =  N1 

CC  4 

>N=N 

K  =  K- 

YT(I 

CCNT 

PETU 

END 


NC  TRANSFORMATIONS  REQUIRED 
LCfc-FASS  TO  LCW-PASS  TRANSFORMATION 
LCU-PASS  TC  HIGH-PASS  TRANSFORMATION 
LCfc-FASS  TO  8ANC-PASS  TRANSFORMATION 
LCV.-FASS  TC  EANC-STCP  TRANSFORMATION 

FFICIENTS  OF  NUMERATOR  POLYNOMIAL 

FFICIENTS  OF  CENCMINATCP  POLYNOMIAL 

ER  OF  NLMERATCR  PCLYNCMAL 

ER  CF  CENOMINATOR  POLYNOMIAL 
TRANSFORMED  NUMERATOR  AND  DENOMINATOR  POLYNOMIALS 

NS1CN  X(25>,Y(25),XT(25),YT(25> 

-M 
+  1 
+  1 
TYFE.NE.l)  GC  TC  11 

1  =  1, M 
1-1 

)  =  >  (I  )=MF1**XM) 
INLE 

1  =  1, M 
1-1 

)  =  t  (I  M(F1**XN  ) 
INLE 
0  2C 
TYPE  .NE.2)  GO  TO  2C 


-i  F1**YF)/(F1**XF) 
+  1 

1  =  1,  M 
1-1 

)=  (X  (K  |*(F1**XM  )  )*FACT 

INLE 

+  1 

1  =  1, M 
1-1 
1 

>  =  Y(K  )*(F1**XN  ) 
INLE 
RN 


C 
C 
C 

r> 

C 
C 

c 


SUSRGUTINE  RCCT(N,RR,RI > 


SL3R0LTINE 
X**N+1=C 
N  :  ORDER 
RF  :ARRAY 
RI  :ARRAY 


TC  FIND  THE  PCCTS  CF  THE  EQUATION 

CF  FCLYNOM IAL 

CONTAINING  REAL  PARTS  OF  CALCLLATED  PCCTS 

CCf^TAINING  IMAGINARY  PARTS  OF  CALCULATED  ROOTS 
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CINENSICh  FF(25),RI(25) 

PI=3. 141592 

XN  =  N 

TETA=FI/XN 

DC  1  1  =  1, rv 

XI =1-1 

ARG=((2.C*XI  )  +  l)*TET4 

RR(I }  =  <IN  UFG) 

RI (I  )=CCS(4RG> 

CCMIME 

RETURN 

ENC 
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APPENDIX  B 


COMPUTER  PROGRAM  LISTING 


COMMCN  h( ICG) ,HM(1C0 ,100) ,X<100) ,V ( 1 00 ) , N, P ,NPULS E 
C 

C  SLBRGUTINE  tc  find  the  time  dcmain  resfcnse 

C  A  DIGITAL  FILTER 

C 

C  N  :  FILTER'S  CRCER  ♦  1 

C  A  :  NUMERATOR  COEFFICIENTS  GF  THE  DIGITAL  FILTER 

C  e  :  CENCNINATCR  COEFFICIENTS   (  WHERE  FIRST  COEFFICIENT 

C  WIuL  BE  NCFNALIZED   ANC  WILL  NCT  EE  ENTERED  TC  THE 

C  PROGRAM 

C 

P»l. 

N  =  24 

NFULSE=2C 
WRITEU  ,1C) 
1C    FORMAT <  '1  •  ) 
CALL    If»FlLS 
CALL    FMR> 
CALL     INFLT 

call  cermet 

WRITEU  ,1C  ) 
CALL  FLCT 
WRITE U  ,1C) 
STOP 
ENC 


SLERCLTINE  HMTRX 
C 

C   SL8R0UTINE  TC  FIND  SYSTEMS  TRANSFER  MATRIX 
C 

COMMCN  h(100) tHM( 100,100 )  , X< 100) ,V (100) ,N,P,NPULSE 

NN=N+1 

NV=NFLLSE+1 

CC  1  1  =  1, NN 

CC  2  J  =  1  »NN 

H  N  ( I  ,  J  )  =  C  .  C 

2  CCNTINLE 
1  CCNTINLE 

CO  2  1*1, Nh 
CG  4  J  =  l  ,NN 
IFU.GT.I  )  GG  TO  4 
L-I-J+l 
HM(I,L  )  =  H(J I 
4  CCNTINLE 

3  CCNTINLE 
WRITEU.2C  ) 

2C    FGRMAT(2X,'    HNATRIX     •///) 
3C    FORMAT  (2X,26F5.3// ) 

CC  5  1*1, NN 

WRITEU  ,2C)  (HM(I ,J)  ,J=1,NN) 
3  CCNTINLE 

RETURN 

END 


SL8RCL1INE  INPUT 
C 
C   SLBROUTINE  TC  FORM  THE  INFLT  FULSE  SEQUENCE 


C 


COMMCN  H(1CQ),HM(100,100>,X(10C),V(100)  ,N,P,NPULSE 
NN=N+1 
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NM=NFUISE+1 

CC  1  I=1,NFLLSE 

>(I)=F 

CCNTIME 

CC  2  I»M»  ,NN 

X(I  }=0.C 

CCNTIME 

kRXTEUvlC) 

FCRMAT  (2X,  •  INPUT  VECTOR 

WRITE  U,2C)  (xm,I  =  l,NN) 

FORMAT(2>  ,26F5.3//) 

RETURN 

END 


//) 


C 

C 

c 


SLBRCLTINE  CCNVGL 

SL3R0UTINE  TC   PERFORM  CONVOLUTION  SUMMATION 

COMMCN  M10J),HM(100,130),X(100),V(100)  ,N,P,NPULSE 
NN=N*1 

co  l  i*i,m 

V(I  )=0.C 
CC  2  J  =  l,Nf* 

V( I  )  =  V(I  )+hMI  ,J)*X(J) 
2  CCNTIME 


1C 
2C 


1  CCNTIME 
WRITE  U  ,1C) 
FORMAT  (2*,'  RESULT  VECTOR 


WPITEU  ,2C  )  <V(I)  ,1  =  1, NN) 
FORMAT  (2>  ,26F5.3//) 
RETURN 
ENC 


//) 


C 

c 

c 


SLERCL7INE    FLOT 

SL8RGUTINE    TC    PLOT    THE    CUTPUT 

COMMGN    H(10O),HM(lOO,100),X(100),V(L0O),N,P ,N PULSE 

CIFEhSKh    T  (250)  ,Y(25C)  ,Z(  25C  ) 

NN=N+1 

NCELTA=1C 

NFCINT=(NFLLSE+3)*NCELTA 

CO     1    I=1,NFCINT 

XI  =  1 

T ( I >  =X  I 

Yd )=C.C 

za  >=cc 

1  CCNTINLE 
K  =  C 

CO    2     I=1,NFCINT,NCELTA 
IF(K.G1.NN  1    GC    TG    5 
K  =  K+1 

CC    3    J«l,4 
L=I+J-1 
Y  (L)=V  (K  ) 
3    CCNTINLE 

2  CCNTIME 

I    CALL     FLCTF(T,Y,NPCINT,1) 


FLCT    INFLT    FLLSES 

K  =  0 

NM=NFLLSE-1 

CO    10     I=1,NPCINT,NCELTA 

[F(K.G1  .hf>  )    GC    TC    7 

K  =  K  +  1 

CC    8    J=l,4 

L=I+J-1 
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Z(U=X  (K  I 

6  COM  I  ME 

1C  CCNTINLE 

7  CALL  FLCTF(T,Z,NPCINT,3) 
RETURN 

ENC 


C 

: 
c 
c 


SL8RCLTINE  IMPULS 

SIBROLTINE  7C  FIND  IMPULSE  RESPGNSE 
(  SOLVING  CIFFERENCE  ECIATIONS  ) 


OF  THE  FILTER. 


B(l4c 

00) 


c 
c 
c 


c 
c 
c 


CGMMCN  H(  ICC) .n 

cimensicn  a(iog> 
cimensicn  clmmy< 
cimensicn  c(10c) 

REAC<5  ,3C  )  N 
NN=N+1 

NNN=NN*2C 


INITIALIZATION  CF  VECTGFS 

CC  35  1  =  1, NNN 
C(I)=C. 
MX)*G« 
MX)*0. 

E(I  )=C. 
35  CGNTINLE 
C(20)=i. 

SUFT  ORIGIN  TC  20 


HM(ICO.IOO) , X(100) ,V(1C0),N,P,NPULS 


3C 
2C 


C  " 


6C 
4C 


U 


NQRDER 

PEAC(5 

REA0(5 

FORMAT 

FORMAT 

WRITEi 

FORMAT 

WRITE ( 

FORMAT 

WRITE( 

wRITE< 

FORMAT 

WRITE ( 

CO  4C 

FACT1= 

NK=II+ 

CO  5C 

NK=NK- 

FACT1= 

CGNTIN 

FACT2= 

NK  =  II 

CO  6C 

M<=NK- 

FACT2= 

CCNTIN 

H(  1 1  )  = 

CONTIN 

h  B  I T  E  ( 

FORMAT 

WRITE* 

FCRMAT 

PQRMAT 

KKK  =  1 

CO  55 

A(KKK  ) 

E(KKK  ) 

KKK  =  KK 


=  M4^ 

,2C  ) 

,2C) 

(12) 

(SF1 

6,1) 

<2X, 

t  ,2) 

(2>, 

6,3  ) 

6,4) 

(2X, 

6,3) 

11  =  2 

C  . 

1 

IL  =  2 

1 

FACT 

LE 

C. 

IL  =  2 

1 

FACT 

LE 

FACT 

LE 

6,21 

(2X  , 

6,  1C 

(2X  , 

(2X  , 

KK  =  2 
=  A  (K 
=  E(K 

K  +  l 


2C 


(Ad  l,I=20,NCRCER) 
(E(I)  ,I=2C,I^CRCEF) 

C.5) 

•  ORDER  OF  THE  FILTER  IS  s  •*I2///I 

•NLMERATCR  CCEFFICIENTS  •///) 
(A(I >,i=2C,NCRCEP) 

'CENOMINATCR  COEFFICIENTS  •///) 

(Ed),  I=2C,NCRCEP) 

CNNN 


CNCRDER 
1+A<IL)*C(NK) 

1,N0RDER 

2*8( IL)*H(NK) 

1-FACT2 

) 

'WEIGHTING  SECUENCE'///) 
KKIlt  I  =  2C,NNN  ) 
12F10.5//) 
12F10.3) 

C  ,NORDEP 
K  ) 
K  ) 
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55  CONTINUE 
K  K  =  1 

CO~65  I=2C,NNN 
CUHMY(KK)=MII 

KK=KK+1 
65  CCNTINIE 

CO  75  l=ltNN 
h(IJ=CLNN\ (I) 

75  cqntinle 
xcum-eo  ) 

B(3) =E(1 > 
E(1)=XCIN 
PEAC(5,1CC  I    F1,F2 
P5A0(5,1CC)  T 
IOC  FGPMAT<2F1C.5> 
NF=1CC 
>NP=1CC. 

FCELTA=(F2-F1)/XNP 

CALL  ZCMN(AfB,M,M,Fl » F2 , FDELTA  ,NP ,T ) 
PETURN 
ENC 
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